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ADIABATIC HEATING OF CONTRACTING TURBULENT FLUIDS 
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ABSTRACT 

Turbulence influences the behavior of many astrophysical systems, frequently by providing non-thermal pres- 
sure support through random bulk motions. Although turbulence is commonly studied in systems with constant 
volume and mean density, turbulent astrophysical gases often expand or contract under the influence of pres- 
sure or gravity. Here, we examine the behavior of turbulence in contracting volumes using idealized models 
of compressed gases. Employing numerical simulations and an analytical model, we identify a simple mech- 
anism by which the turbulent motions of contracting gases "adiabatically heat", experiencing an increase in 
their random bulk velocities until the largest eddies in the gas circulate over a "Hubble" time of the contraction. 
Adiabatic heating provides a mechanism for sustaining turbulence in gases where no large-scale driving exists. 
We describe this mechanism in detail and discuss some potential applications to turbulence in astrophysical 
settings. 

Subject headings: hydrodynamics — turbulence 



1. INTRODUCTION 

Turbulence - the bulk random motion of a gas or 
fluid - is ubiquitous in astrophysics. Turbulenc e can 
be ge nerated by instabilities, including gravitational (iJeansI 
[19021) . s hear (vonHelmholtz 1868; K elvinl ITSTlh . con- 
vective (jRavleigh 1884; Tavlon il95Q|) . and magnetoro- 
tational i Balbus & Hawley 199 T|). For o verviews, see 
lElmegreen & Scalo 



(2004) and McKeE&Ostriker (2007). 
Given its wide-spread importance, turbulence remains a criti- 
cal area for astrophysical research. 

Supersonic turbulence has been studied in great de- 
tail with numerical simulations. For instance, super- 
sonic isothermal turbulence exhibits a lognormal density 
distribution, with a w idth that increases with the Mach 
number (e.g., IVazquez- Semadeni 1994; Padoanetal. 1997c 
iKritsuk et al.l 120071: iLemaster & Ston e 2008; F ederrath et al J 
|2010; Price et al. 20TTh. The properties of supersonic isother- 
mal turbulence appear to be independent of the s imulation 
methodology (e.g., Kitsionas et al. 2009: [Pnce & Federrath 
HpiO; Bauer & Springel 2011). Most studies of turbulence 
have involved gases simulated in a static volume, whereas as- 
trophysical gases often expand or contract under the influence 
of pressure or gravity. Little is currently known about the de- 
tailed structure of expanding or contracting turbulent gases. 

In this Letter, we examine the behavior of turbulence during 
the contraction of a gas arising from pressure or self-gravity. 
In Section [21 we use simulations to model contracting turbu- 
lence and demonstrate that turbulence adiabatically heats dur- 
ing contraction provided the eddy turnover time[j is shorter 
than the contraction time. We term this mechanism "adiabatic 
heating", and in Section[3]we present an analytical model that 
successfully describes its behavior. We discuss some poten- 
tial astrophysical applications of adiabatic heating in Section 
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^ While the term eddy accurately describes the vortices of incompressible 
turbulence, it is less accurate for motions in compressible turbulence. Lack- 
ing a better term, we nonetheless refer to the large scale motions in compress- 
ible turbulence as eddies. Similarly, the term turnover time is used to describe 
the timescale of these motions. 



[H and summarize and conclude in Section [51 
2. NUMERICAL SIMULATIONS OF CONTRACTING TURBULENCE 

We use hydrodynamic simulations to study turbulence in a 
contracting background. We model the contraction by param- 
eterizing the changing physical size l(t) and coordinate scale 
factor a = lit) /L of di cubic volume of initial length L through 
a "Hubble" parameter H = d/a that may depend on time t. In 
terms of the proper coordinates within an isotropically con- 
tracting volume, the Euler equations connecting derivatives 
of the density p, momentum pv, and pressure p are altered by 
terms that depend on the Hubble parameter as 



dpv 



^=-S/-(pv)-3Hp 
= -V • (pvv) -Vp- AHpv 



(1) 



(2) 



(see, e.g.. Section 9 of lPeeblesll l98Q'). For an initially constant 
density and velocity gas without dissipation, these terms give 
rise to two important scalings well known from cosmology: 
p (X oT^ and v (xaT^. Mass conservation dictates the density 
scaling, but the presence of dissipation (either through phys- 
ical viscosity or numerically through the discretized form of 
Equation [2]) implies that the adiabatic velocity scaling does 
not strictly hold in the contraction of a turbulent gas. How 
the turbulent velocity evolves depends on how dissipation op- 
erates during the contraction, and simulations are required to 
provide a detailed description. 

The simulations were performed using a version of the 
magnetohydrodynamicscodeA^/z^^a (Stone et al. 2008) mod- 
ified to model contracting and expanding turbulent gases (see 
Equati ons [T] and [2l and below). Athena is a grid code based 
on the lGodunovl ^[959) meth od. The calculation s use piece- 
wise parabolic reconstruction (IColella & Woodw ard 1984) to 
extrapolate initial states for the Riemann problem between 
cells and compute final states using an exact solver (Torol 
Il999|) . Cell-aver aged conserved quantitie s are updated using 
unsplit methods (|Gardiner & Stond [2008l). Our modifications 
to Athena include a Runge-Kutta integrator to evolve a dif- 
ferential equation for the scale factor a that depends on the 
possibly time-dependent Hubble parameter H. 
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Figure 1. Contraction of a turbulent gas and the adiabatic heating mechanism. Shown are the RMS turbulent velocities for contracting isothermal simulations 
with differing Hubble parameters (colored lines). In an initially "slow" contraction (top left panel) where the contraction frequency |//| is lower than the eddy 
turnover frequency cj, large-scale eddies circulate and the turbulent cascade transfers energy to small scales where it dissipates. Instead, if the contraction 
frequency is initially higher than the eddy turnover frequency (top center and right panels), the turbulence heats roughly adiabatically at first. As the ratio of 
eddy turnover and contraction frequencies in the simulations become comparable, whether the turbulence heats or decays as the scale factor « — )► depends on 
the evolution ofH{a). In a dynamical contraction (top center panel) where H oc a~^/^ the turbulence continues to heat after uj tracks \H\, whereas for constant H 
(top right panel) it decays. The solid black lines in each panel show the evolution of turbulent velocities predicted by our analytic adiabatic heating model (see 
Section[3). The dotted lines indicate the model predictions for the synchronized ratio of u;/\H\. To indicate the effect of the forced turbulence initial conditions 
on the contracting turbulence, the Simulation A panels show two realizations (thick and thin lines). 



The initial conditions are snapshots of driven isothermal 
turbulence (with sound speed c^ = 1 and mean density p=l) 
simulated on an N = 512^ re solution grid. The r andom forcing 
field is generated following Bertschinger ("2001), with power 
input into the two larges t modes in the unit (L = I) peri- 
odic box (e.g.. iKritsuk et al.ll200 7). Driving at intermediate 
scales produces similar results. A Helmholtz decomposition 
in Fourier space removes the dilatational component and each 
forcing field is normalized to maintain an average Mach num- 
ber M ^ 6 when applied as an acceleration ten times per cross- 
ing time Across ^ L/2Mcs. The driving is applied for ten cross- 
ing times and then terminated before the contraction initiates. 
Although we simulate isothermal gases, results relevant for 
the adiabatic heating mechanism originate from Equations \T\ 
and[2]and should generalize to other adiabatic indices. 

2.1. Models of Contraction 

Two rates characterize the isotropic contraction of a turbu- 
lent gas, the contraction frequency \H\ = \d/a\, also called 
the Hubble parameter (// < for a contraction), and the eddy 
turnover frequency 

uj^—-, (3) 

aL 

where v is the root-mean- squared (RMS) turbulent velocity 
(for isothermal turbulence v = Mc^, where M is the typical 
Mach number and Cs is the sound speed). In our simulations, 
the values of v(^ = 0) ^ 6c^, c^ = 1, and L=l imply an initial 
eddy turnover frequency uj(t = 0) ^ 6. To demonstrate the 



generality of the adiabatic heating mechanism we simulate 
three different scenarios for the time-dependent relation 
between// andcj: 

Simulation A: an initially "slow exponential contrac- 
tion", with H = constant. In this case the scale factor evolves 
as a(t) = aoQxp[H(t-to)], where ao = 1 is the initial scale 
factor and to is the time when the contraction ensues. We 
choose a constant Hubble parameter H = -2, such that the 
contraction is slow (i.e., \H\ <C cj) initially. To characterize 
run-to-run variations, we perform two such simulations 
differing only in their forced turbulence initial conditions. 

Simulation B: an initially "fast dynamical contraction", 

with H (X ^/p. In this case the contraction time ^cont ~ I^T^ 
scales with the dynamical time ^dyn oc l/v^ set by the mean 
density p (xa~^. Starting with an initial value H(t = to) = Ho, 
the Hubble parameter varies with the scale factor as 
H = Ho(a/ao)~^^^. The scale factor decreases with time as 
a(t) = ao[3Ho(t-to)/2-\-l]^/^ (recall that Hq < 0). We set 
Hq = -20 to induce an initially fast (|//| ^^ uj) contraction. 

Simulation C: an initially "fast exponential contrac- 
tion", with H = constant. In this case the scale factor evolves 
with the same time-dependence as in Simulation A, but at a 
constant contraction frequency (\H\ = 20) such that \H\':^ uj 
initially. 
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Figure 2. Adiabatic heating of isothermal supersonic turbulence during contraction. Shown is the logarithmic density distribution through a thin (8/512) slice of 
Simulation B at three values of the scale factor a. Here, the Hubble parameter scales with the inverse dynamical time of the gas. The density maximum in the 
color map of each panel is scaled by a~^, and the dynamic range of each image is 10^. The turbulence adiabatically heats from the initial Mach number of M ~ 6 
to M ~ 1 1 as the scale factor decreases by a factor of ~ 2. The presence of limited dissipation breaks the perfect adiabatic scaling. The bulk properties of the gas 
in each panel, including the density distribution and intermittency, behave similarly to static-frame isothermal turbulent gases at the same Mach number. 



2.2. Simulation Results 

The properties of contracting turbulence evolve with de- 
creasing scale factor a in a manner that depends upon the ratio 
of the eddy turnover frequency to the contraction frequency 
(Figure [B. The simulations demonstrate that if the contrac- 
tion is slow (\H\ <^uj initially in Simulation A, left panels), 
the turbulent velocities decay, whereas if it is fast (\H\ ^ uo 
initially in Simulations B and C, center and right panels), the 
turbulent velocities amplify. In a slow contraction, large vor- 
tices circulate and nonlinear interactions transfer energy to 
smaller scales where it is dissipated before the box shrinks 
appreciably. When the contraction is fast, energy bearing ed- 
dies are adiabatically compressed, dissipation primarily op- 
erates on small scales, and turbulent velocities increase. In 
each example in Figure[T] uj/\H\ become comparable (bottom 
panels). From this trend we surmise that the eddy turnover 
frequency may eventually "synchronize" with the contraction 
frequency. Our simulations provide a hint of this behavior, but 
the synchronized state is not well-explored. Physical consid- 
erations suggest the synchronization is stable since the eddies 
are compressed on their circulation timescale, and the veloc- 
ities of large eddies should hover around v ~ \H\aL. If true, 
as a ^ 0, the velocities would decrease for constant \H\ as 
observed in Simulations A and C but continue to increase for 
\H\ (X a~^/^ as seen in Simulation B. 

Figure [2] shows the projected density distribution through 
a slice of Simulation B at three scale factors during the con- 
traction. Initially, at a = 1 (left panel) the gas has a turbulent 
velocity v = Mcs ~ 6c^, a lognormal density distribution, and 
a velocity power spectrum characteristic of supersonic turbu- 
lence. Compression during contraction heats the turbulence 
in the isothermal gas to Mach number M = 1 A (M = 11.2) 
by scale factor <2 = 0.81 (a = 0.43). As the Mach number in- 
creases, the width of the lognormal density distribution in- 
creases, the intermittency amplifies, and the velocity power 
spectrum steepens much in the same way as driven, constant 
volume isothermal turbulen ce simulations be have as a func- 
tion of Mach number (e.g., iPrice et al.ll2011h . The increase 
in turbulent velocities arises from the approximate inverse de- 
pendence of velocity on the scale factor during contraction. 
Since a turbulent cascade transfers large-scale power to small 
scales where it dissipates, the degree to which the turbulent 
velocity tracks a~^ during the simulation depends on the rate 
of energy transfer to small scales. 



3. ANALYTICAL MODEL OF THE ADIABATIC HEATING 
MECHANISM 

We model the behavior of turbulent velocities during con- 
traction by calculating the approximate time rate of change 
of the kinetic energy per unit mass in the gas, including two 
important terms. A term capturing the effects of adiabatic 
heating follows from noting that the adiabatic velocity scal- 
ing implies va = vq, where vq is a constant. Thus 



d 



AH 



dv voda 2 rA^ 

= y—- = -v-r- -- = -if v^. (4) 

dt a^ dt 



A second term capturing the rate of kinetic energy dissipation 
is modeled using a parameter r] that describes the efficiency 
of the energy cascade. Physically this term would represent 
viscosity in the Navier-Stokes equation, but in our simulations 
it arises from dissipative truncation error in the discretization 
of Euler's equations. The dissipative term reads 



(5) 



<^t\2j,^,, \l' 



The relevant length scale in Equation [5] is the driving scale 
(e.g.. iMac Lowl [19991) , which is the box size a L in our cal- 
culations. Simulations of driven incom pressible (jGotoh et al.l 
I2002t lBeresnvakll2QTTh and transonic (iSchmidt et al.l I2006h 
turbulence suggest that r] ^ I. The total rate of change in 
the turbulent velocity is then 



dv V 

= -Hy-r]—. 
dt aL 



(6) 



Noting that aH = da/dt, and writing the eddy turnover fre- 
quency as uj(Y^a) ~ y/aL, the rate of change of the turbulent 
velocity with scale factor can be recast as 



dv / uj\y 
da~~\ '^^HJa' 



(7) 



We will refer to Equation |7] as the "adiabatic heating equa- 
tion" and it provides quantitative insight into the slow and fast 
contraction regimes discussed qualitatively above. It shows 
that the heating of turbulence during contraction is moderated 
by dissipation with an efficiency proportional to ^ uj/\H\. 
When the contraction is slow, the dissipative term is larger 
than the heating term (r]uj/\H\ :>l) and the velocity decreases 
with decreasing scale factor (dy/da > 0). Conversely, when 
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the contraction is fast, the velocity increases with decreasing 
scale factor (dv/da < 0). Once the eddy turnover and collapse 
frequencies become comparable and synchronize, whether the 
velocities grow or decay depends on how H varies with a. 

Figure [T] depicts the evolution of the turbulent velocity v 
(upper panels) and the ratio of frequencies uj/\H\ (lower pan- 
els) obtained from simulations, along with predictions from 
the adiabatic heating equation using a dissipation parameterQ 
fixed at 7/ = 1.2. The model reproduces well the behavior of 
both properties of turbulence in contracting gases. For simula- 
tions where oo/\H\ ^ 1 ox uo/\H\ <C 1 during the whole com- 
putation (left and center panels) the adiabatic heating equa- 
tion provides an accurate description of the turbulent velocity 
evolution. In the initially fast contraction with constant H 
(right panel) the general behavior is also well modeled by the 
adiabatic heating equation, but the predicted transition from 
heating to dissipation occurs later in the model than in the 
numerical computation. These differences may indicate that 
dissipative effects are delayed in the simulation by an eddy 
turnover time relative to the analytical model. 

The model suggests the adiabatic heating mechanism drives 
the turbulence in the contraction to an asymptotic relation be- 
tween uj and \H\ determined by the scale factor dependence 
of the Hubble parameter. From the adiabatic heating equation 
and the definition of the eddy turnover frequency, we obtain 



dlog(uj/H) 
dlog(l/a) 






1 



dlogH 
dlog(l/a)' 



(8) 



The asymptotic relation is approached as 
dlog(uj/H)/dlog(l/a) -^ 0. For H = constant we ex- 
pect cjo/\H\ -^ 2/r] (Simulations A and C), while for 
H oc a-^/2 we have uj/\H\ -^ l/2r] (Simulation B). We find 
that the simulations follow the evolution in uj/\H\ predicted 
by Equation [U (see Figure 1, lower panels), but can show 
substantial run-to-run variations (e.g.. Simulation A). 

The asymptotic relation between the typical turbulent ve- 
locities and scale factor can be deduced from Equation [8] By 
defining 



dioga 



(9) 



from the adiabatic heating equation we find simply that v ex 
a^~^ once the contraction and eddy turnover frequencies have 
synchronized. 

Although the turbulent velocity roughly tracks the expected 
scaling v oc a~^ before an eddy turnover time elapses, the de- 
gree of adiabaticity depends on the small scale dissipation 
rate. In additional simulations of contracting incompress- 
ible isothermal turbulence (M ~ 0.05), we have found that 
even low resolution simulations display almost exact adia- 
batic heating. For very large isothermal turbulent velocities 
(M > 10), we have found that the heating becomes more adi- 
abatic with increasing resolution. This sensible behavior does 
not affect any of the presented results which have been tested 
over a wide range of resolutions (from N=64^ to N=512^) and 
display similar behavior for the same choice of Hubble pa- 
rameter evolution and initial turbulent velocity. 

4. DISCUSSION 

Our results have many potential astrophysical applications, 
but bear especially on the problem of turbulence in giant 

^ We find that using ry = 1.2 in Equation |7]reproduces well the isothermal 
simulation results at all Mach numbers. 



molecular clouds (GMCs). The properties of GMCs are 
likely set by turbulence, as the re lations betwe en cloud ve- 
locity dispersion, size, and mass (lLarson|[T981h may reflect 
properties of tur bulence through the vel ocity structure func- 
tion (see, e.g., Elmegreen & Scalo 2004). If so, observations 
of the dispersion- size re lation of GMCs (e.g., Bolatto et al.l 
120081: iHeyer et aH |2009|) ro ughly agree with the properties 
of compressible turbulence dB allester os-Paredes et al.l 120061: 
iKritsuketaP 120071: iFederrath et al. 20l3). It has been ar- 
gued that molecular clouds are in approximate virial equi- 
librium (Larson 1981; Solomon etal. 1987) such that tur- 
bulent motions either balance the cloud self-gravity or oth- 
erwise reflect the depth of the gravitational potential (e.g., 
iBertoldi & McKedl 19921) . This apparent virial balance poses 
a significant challenge for understanding the origin and 
evolution of turbulence in molecular clou ds, as turbulence 
should dissipate on a crossing time (e.g. JGoldreich & KwanI 
1974) that is short er than estimates of the cloud lifetime 
(iBlitz & Shul 119801) . The typical GMC Hfetime is de- 
bated because they may not be in exact balance or grav- 
itationaUy bound (iHartmann et al.l 120011: iDib et al.1 120071: 
Dobbs et al. '2011), and perhaps undergo frequent collisions 
(e.g., Tasker & Tan 2009). It remains unclear how turbulence 
could generically provide support against gravit ational col- 
lapse since without driving it quickly dissipates (IStone et al.l 



119981: Mac Low etal.ll 19981: ICho & Lazarianll2003l) 

Our study suggests that the connection between velocity 
dispersion and size may reflect the competition between adi- 
abatic heating and dissipation. Depending on the nature of 
the contraction, the adiabatic heating mechanism can enable 
the typical turbulent velocity to scale with a positive power 
of the size of a contracting cloud or region, preserving a 
connection between velocity dispersion and cloud size with- 
out an external source for driving the turbulence. From the 
discussion in Section [3l the observed scalings of v ex L^/^ 
(ISolomonetal.l[T987l: iHever et al.l l2QQ9l) require H ex a'^^^. 
This scaling is quite different than what might occur in a 
gravitational collapse, where naively one expects H ex a~^^^ 
and V ex L~^/^. If turbulent velocities in GMCs do not orig- 
inate from gravitational collapse, adiabatic heating may still 
provide a method for instilling the observed scaling relations 
through other compression mechanisms. 

Although our study has focussed on adiabatic heating in 
isotropically contracting turbulence, we have examined other 
scenarios. The inverse process ("adiabatic cooling") similarly 
operates in expanding gases. Using simulations with a Hubble 
parameter if > 0, we have verified that expanding turbulent 
gases adiabatically cool if the eddy turnover frequency is less 
than the expansion frequency H. If the expansion is very rapid 
(H :^ uj) then the turbulence freezes out with v ex a~^ Adia- 
batic cooling may be relevant for turbulent astrophysical sys- 
tems that rapidly expand, such as those formed in high speed 
impacts or explosions. We have also simulated anisotropic 
systems, and found that adiabatic heating and cooling can op- 
erate simultaneously in different directions depending on the 
sign of the effective Hubble parameter for each axis. Appli- 
cations of anisotropic comp ressions include studies of shock- 
turbulence interaction (e.g. jAdams & Shariffll 19961) . 

Some previous works presented ideas related to adiabatic 
heating. Olson & Sachs (1973) studied analytically the evo- 
lution of mean vorticity in incompressible turbulence in ex- 
panding universes without dissipation, and commented that in 
a contracting universe the vorticity would "blow up". There 
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were early works e xploring whether turbulenc e could seed 
structure formation (lJoneslll976t lQzemoiiri978|) that consid- 
ered the adiabatic scaling of velocity with inverse scale factor 
The collisional N-body calculations of IScalo & Pumphreyl 
([19821) suggested that turbu lence m ight slow a gravitational 
collapse. Vazquez-Semadeni eLaU {1998) suggested that tur- 
bulent velocites might depend on the mean density during col- 
lapse. Our description of adiabatic heating has combined and 
expanded upon some of these concepts. 

5. SUMMARY 

Using simulations of contractin g isothermal turbu lent gases 
performed with the Athena code dStone et al.ll2008|) . we have 
identified an "adiabatic heating" mechanism by which ran- 
dom bulk motions are amplified by compression. Adiabatic 
heating acts to increase the turbulent velocities of a contract- 
ing gas if the frequency (or Hubble parameter) \H\ of the con- 
traction is larger than the eddy turnover frequency u ~ w/aL 
(L is the initial box size, a is the scale factor of the contraction, 
and V is the turbulent velocity). When \H\:^uj the cascade of 
energy from large scales to small scales is limited and dissipa- 
tion becomes inefficient, thereby allowing the gas velocities to 
heat with the adiabatic scaling y (xa~^ expected from Euler's 
equations in a contracting background. When \H\ <^uj, the 
cascade operates efficiently and energy dissipation proceeds 
similarly to turbulent decay in static volumes. In each case, 
the turbulent velocities evolve toward cj/\H\ ~ 1 during con- 
traction. Using these insights, we develop an analytical model 
to describe the rate of change of energy per unit mass in the 
gas as a competition between compressive adiabatic heating 
and dissipation on small scales. The analytical model suc- 
cessfully predicts the dependence of both the RMS turbulent 
velocity v and the frequency ratio cc;/|^| on the scale factor a. 

BER is grateful for generous support from the University 
of Arizona and Steward Observatory. PMG thanks Jungyeon 
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cluster at the California Institute of Technology Division of 
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